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The seawater Mg/Ca ratio increased significantly from ~ 80 Ma to present, as suggested by studies of 
carbonate veins in oceanic basalts and of fluid inclusions in halite. We show here that reactions of 
mantle-derived peridotites with seawater along slow spreading mid-ocean ridges contributed to the 
post-Cretaceous Mg/Ca increase. These reactions can release to modern seawater up to 20% of the yearly Mg 
river input. However, no significant peridotite-seawater interaction and Mg-release to the ocean occur in 
fast spreading, East Pacific Rise-type ridges. The Mesozoic Pangean superocean implies a hot fast spreading 
ridge system. This prevented peridotite-seawater interaction and Mg release to the Mesozoic ocean, but 
favored hydrothermal Mg capture and Ca release by the basaltic crust, resulting in a low seawater Mg/Ca 
ratio. Continent dispersal and development of slow spreading ridges allowed Mg release to the ocean by 
peridotite-seawater reactions, contributing to the increase of the Mg/Ca ratio of post-Mesozoic seawater. 



Has the composition of seawater changed through time? It appears that the Mg/Ca ratio of the Mesozoic 
ocean was 3 to 5 times lower than that of modern oceans 110 (Fig. la). Seawater Mg/Ca ratio is important as 
it affects the growth of calcitic versus aragonitic marine organisms, the deposition of inorganic carbonates 
(ooids, cements), as well as ocean-atmosphere C0 2 exchange, which influences climate. 

A number of explanations have been suggested for the increase of the Mg/Ca ratio in post-Mesozoic seawater. 
One calls for a post-Mesozoic decrease of both the rate of Mg removal and of Ca release by the basaltic crust due to 
an alleged decrease in the rate of seafloor spreading and of hydrothermal circulation 1112 . Another calls for 
enhanced deposition of dolomite that contributed to the late Mesozoic low seawater Mg/Ca ratio, followed by 
a decrease in dolomitization during the last —100 Ma, due to lowering of sea level that restricted shallow seas 
where dolomitization is favored 1,6,7,9 . Post Mesozoic increased continental weathering has also been suggested 13 , as 
it would increase the river input of Ca and Mg followed by preferential deposition of Ca as carbonates resulting in 
a Mg/Ca increase. These processes, although they can contribute to solve the problem, fail to explain satisfactorily 
the strong post-Mesozoic increase of seawater Mg/Ca 9,10 . We suggest here that seawater-mantle peridotite reac- 
tions in slow-spreading ridges may provide an important additional mechanism. 

A number of processes affect the Ca and Mg content of the oceans. Present day average river input derived from 
continental weathering has been estimated at ~ 13.2 - 10 12 to ~ 15 - 10 12 mol/yr Ca and at ~ 5.2 - 10 12 to ~ 
6.1 • 10 12 mol/yr Mg 1,6,14 . Before the discovery of subridge hydrothermal circulation no mechanism was known 
to get rid of Mg brought into the ocean by rivers 14 . During subridge high-T seawater circulation in the basaltic 
crust some Ca is extracted from basalt and added to seawater; in contrast, Mg is lost by seawater and incorporated 
into secondary minerals in the crust 15 . The quantity of seawater that goes through high-T hydrothermal circula- 
tion within the basaltic crust along today's mid-ocean ridges has been estimated at ~ 3-5 - 10 13 kg/yr (ref. 16). 
Geochemical (Li and Tl isotopes) and geophysical data 17,18 limit the high-T hydrothermal flux from the oceanic 
crust at ~ 1-3* 10 13 kg/yr, in contrast with a flux of ~ 5-10* 10 13 kg/yr necessary to balance the Mg oceanic 
budget 19 . Off- axis low-T alteration of the basaltic crust may affect significantly the budget of oceanic elements 20-22 , 
including Mg, although Tl geochemistry 18 and estimates of off-axis basalt alteration 23 suggest little low-T basalt/ 
seawater chemical exchange at ridge flanks. Uncertainties in off-axis water fluxes 18,20 and variability in the 
chemistry of low-T off-axis hydrothermal fluids 21,22 make it difficult to estimate Mg fluxes at ridge flanks. 
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Figure 1 | Seawater concentration of Mg and Ca. (a), Estimates of Mg/Ca ratio of seawater since 180 Ma from ref. 10, including data from halite fluid 
inclusions (circles), ridge flank basalt carbonate veins 10,71 (triangles), echinoderm ossicles (green diamonds), pore-fluid chemical profile modelling 
(yellow diamonds), benthic and planktonic foraminifera (grey and white squares), and rudists 72 (yellow squares). The Mg/Ca temporal variations 
(red circles and diamonds) are obtained from models of secular variations of seawater concentration of Ca and Mg shown in (b) and (c), respectively, 
(b) , Changes through time of seawater concentration of Ca versus variations in spreading rates and accretionary boundary geometry over the last 150 Ma, 
assuming an initial concentration of 24 mmol/kg (from brine inclusions in 150 Ma old evaporites 5 ), and assuming constant river influx and residual 
outflow; variable Ca-release flux from hydrothermal circulation in basalt; and different degrees of Ca-removal variable-flux due to alteration of mantle 
rocks exposed at the seafloor (red filled circles and open diamonds), (c), Models of Cretaceous and Cenozoic changes in seawater concentration of Mg due 
to variations through time of oceanic crustal production and volume of mantle that can react at low-T with seawater (50%, red circles, 0% and 100% red 
diamonds), assuming constant river influx; variable Mg-release flux from low-T peridotite-seawater reactions; variable Mg-removal flux from high-T 
hydrothermal circulation in basalt; and a constant residual outflow, starting from an initial concentration of 30.5 mmol/kg at 150 Ma 5 and assuming a 
total seawater mass of 1.33 8 1 0 21 kg. Seawater Mg concentration due solely to Mg-release from MORP-seawater reactions for different volume fractions of 
peridotite interacting at low-T with seawater (25%, dashed-dotted line; 50%, solid line and 100% dashed line) is also shown. Red square is the modern 
seawater concentration. Blue and cyan circles are ancient seawater compositions calculated from measurement of fluid inclusions in marine 
evaporites 3,5,73 . (d), Triangular Mg-Si-Fe diagram showing the bulk composition of hydrated mantle-derived peridotites (open symbols) and 
reconstructed primary compositions of the same rocks (filled symbols). Diamonds refer to SWIR peridotites 34 , circles to Vema Lithosperic Section 
peridotites 38,39 . Both populations show Mg-loss in serpentinites compared to unaltered primary peridotites. 
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Given a concentration of Mg in seawater of 0.053 mol/kg, and 
assuming high-T hydrothermal fluxes of ref. 16 extracting Mg totally 
from seawater, it follows that ~ 1.6-2.7 - 10 12 mol/yr Mg are lost by 
seawater due to high-T hydrothermal flow. This Mg loss is a signifi- 
cant fraction of the yearly Mg river inflow into the oceans 20,24 . In 
addition, high-T hydrothermal reactions transfer ~ 0.8-1.25" 10 12 
mol/yr Ca 6,15,20,24 from basalts to seawater. 

Precipitation of carbonates from sea water is affected by temper- 
ature and seawater Mg/Ca ratio 2,7 , with high Mg/Ca ratios favoring 
aragonite over calcite 11 . Deposition of dolomite extracts Mg from 
seawater at the rate of 1.7' 10 12 mol/yr according to ref. 6. 

Results 

Our suggestion that mantle peridotite-seawater relations have con- 
tributed to the post-Mesozoic increase of seawater Mg/Ca ratio is 
based on the following steps (see Supplementary Information for 
detailed explanations): (i) we modeled mantle dynamics beneath 
the global ridge system, showing that mid-ocean ridge peridotite 
(MORP) distribution depends on spreading rate. We assessed the 
volume of MORP that has the potential to react with seawater under 
conditions allowing Mg release; (ii) we demonstrated that MORP 
samples have lost Mg relative to their primary unaltered parent; 
(iii) we showed through global plate reconstructions how crustal 
production and average spreading rate have changed since the 
Cretaceous; (iv) we combined the results of (i), (ii) and (iii) to cal- 
culate the fluxes of Mg and Ca to the oceans as a result of mantle 
rock-seawater interactions, (v) we combined these results with esti- 
mates of the riverine and high-T hydrothermal Mg and Ca fluxes to 
model seawater Mg, Ca and Mg/Ca since 150 Ma (Fig. 1). 

Significant stretches of today's slow- spreading ridges have mantle- 
derived peridotites emplaced either on, or at shallow levels below the 
sea floor, so that seawater can have access to the peridotites (Figs. 2 
and 3). Peridotite-H 2 0 reactions have been discussed in a large body 
of literature 25 29 . They can take place at temperatures up to ~ 500°C, 
with or without volume increase and elemental exchange 27 , other 
than acquisition of H 2 0 by the rock. Some reactions allow Mg to 
be released in solution 30 . In contrast to basalt-seawater hydrothermal 
reactions, when sea water reacts with peridotite Mg can be extracted 
from the rock, provided the temperature of the system is < 150°C 
and water/rock (W/R) ratio is high 30 34 . 

Experiments 35 and numerical simulations 36 on olivine-H 2 0 reac- 
tions show that Mg can be released in solution even at 300°C. Mg is 
also released by low-T dissolution of brucite, a phase likely to form 
during higher-T serpentinization, but hardly ever found in ocean 
floor serpentinites 25,30,34 . Low-T incongruent dissolution of olivine 
and enstatite may also release Mg 34,37 . 

Modern MORP contribution to the ocean Mg and Ca budget. In 

order to test if Mg is lost during reactions of seawater with MORP, we 
compared the chemical composition of primary unaltered MORP 
with that of hydrated MORP samples obtained at over 20 sites 
where mantle is exposed along the Vema Lithospheric Section in 
the central Atlantic 38,39 . They represent mantle residual after 
extraction of basaltic melt at a single segment of the Mid Atlantic 
Ridge (MAR) throughout 26 million years of crustal accretion. These 
rocks all contain > 10% H 2 0. Their whole-rock primary major 
element content was calculated by combining their reconstructed 
primary modal composition with the chemical composition of the 
primary phases olivine, orthopyroxene, clinopyroxene and spinel 
(see Supplementary Information). We found that the altered 
peridotites, other than having gained H 2 0 and lost some Si, have 
lost about 5-6% Mg relative to their unaltered parents (Fig. Id). A 
similar loss of Mg due to reactions with seawater has been 
documented also in peridotites from the SW Indian Ridge 
(SWIR) 34 and from the Atlantis Massif (MAR at 30° N) 40 . A 
comprehensive major and trace element study of ~ 130 abyssal 



peridotite samples from Pacific and Indian ocean ridge-transform 
systems shows a mean ~ 10% wt% MgO loss relative to the MgO 
content of unserpentinized protoliths 41 . 

The extent to which ocean floor serpentinization implies volume 
increase is not settled yet. Textural analysis of thin sections of the 




Figure 2 | Mantle seawater reactions at mid ocean ridges, (a), Present-day 
half spreading rates along mid ocean ridges, using relative plate kinematic 
models from ref. 45. Thick red solid lines indicate ridge segments, mostly 
with half spreading rate < 10 mm/yr, that can contribute to <150°C 
MORP-seawater reactions and consequently to Mg release. The global 
ridge system was divided in 80 sectors each 2048 X 1024 km, partially 
overlapped to avoid edge effects. Assuming a plate spreading velocity for 
each sector by averaging spreading rates and directions, we calculated 
passive mantle flow, thermal structure and melt production beneath each 
ridge segment following methods of ref. 44. We calculated along-axis 
crustal thickness assuming pure-fractional melting and complete melt 
extraction. We neglected latent heat of fusion by melt freezing and 
hydrothermal cooling. The along-axis depth-distribution of the 150°C 
isotherm was obtained by averaging depths from 10 km-wide across-axis 
sections. The integral of differences between isotherm depth and crustal 
thickness (where they are positive) times the full spreading rate gives us an 
estimate for each ridge segment of the volume of mantle-derived MORP 
that can interact yearly with seawater at a temperature below 150°C. 
(b), Half spreading rates along mid-ocean ridges at 83.5 Ma, obtained 
using digital models for spreading velocities from refs. 47 and 48. Half 
spreading rate was generally faster than 15 mm/yr, implying little or no 
MORP-seawater reactions. Reconstructed positions of continents at 
83.5 Ma (Santonian), relative to Africa assumed to be fixed, were obtained 
with the Opiates software (http://www.gplates.org). 
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Figure 3 | Effects of spreading rate on the possibility of MORP-seawater reactions during hydrothermal circulation at mid ocean ridges. Assuming a 
constant mantle potential temperature of 1350 °C at a depth of 150 km and plate thickening passive flow, the degree of melting of the upwelling mantle 
depends on whether plate separation is slow (a) or fast (b). We solved for the steady-state three-dimensional passive mantle flow via a Fourier pseudo- 
spectral technique 44 . The base of rigid plates, assumed to correspond to the depth of the 800°C isotherm, was obtained iteratively solving each time the 
mantle temperature field, starting from a constant-thickness plate-flow model. Streamlines show a schematic passive corner flow in the asthenosphere. 
The shaded triangles indicate the fraction of melt generated across axis, including the effect of water on peridotite solidus. Assuming a ~ 0.02 wt% H 2 0 
content in the upper mantle, the peridotite solidus is lowered causing partial melting in a subridge mantle region wider and deeper than if the mantle were 
dry. The faster the mantle rises, the shallower melting ends. As a result, more mantle melts, creating a thicker basaltic crust and shallow seawater 
penetration, with little chance for mantle ultramafics to interact with seawater and to release Mg to the ocean. Slow ridge MORP-seawater reactions may 
result in Mg release to the ocean. 



Vema serpentinites suggest that volume increase during serpentini- 
zation was limited to <20%. Similar low volume increase was esti- 
mated for Atlantis Massif serpentinites 42 , implying significant Mg 
loss 25,30 . Nevertheless, Mg-bearing fluids derived from low-T hydro- 
thermal circulation in mantle peridotites have not been reported to 
date, consistent with bulk- rock Mg depletion being probably related 
to low-T "pervasive weathering" by seawater of a relatively thick 
ultramafic zone 34,41 . However, the issue of Mg-loss versus temper- 
ature during MORP-seawater reactions is constrained by the tem- 
perature below which the reacting fluid is undersatured in Mg-rich 
minerals (i.e., < 150°C following ref. 34), independently of the pro- 
cess leading to Mg-depletion of MORPs (i.e., near constant-volume 
mantle hydration or pervasive "marine weathering"). 

We estimated next the potential for peridotite-seawater reactions 
in the modern versus Cretaceous-Cenozoic oceans (Fig. 2). Seawater- 
MORP reactions require: (a) thin (or absent) basaltic crust; (b) sub- 
seafloor seawater penetration (unless peridotites are exposed directly 
on the seafloor); (c) temperature below the 500°C isotherm ( 1 50°C to 
favor Mg release). 

Basaltic crust thickness is related to subridge thermal structure 
and increases with spreading rate 43 . We calculated it (see Supple- 
mentary Information) following computational methods on plate- 
driven subridge mantle flow for a lithosphere that thickens with age 44 
(Fig. 4). We considered the global mid-ocean ridge system using the 
current plate motion model of ref. 45. Half spreading rates, ranging 
from 3.9 mm/yr (i.e. Gakkel Ridge) to 75 mm/yr (i.e. East Pacific 
Rise), were used to calculate mantle flow velocity. Crustal thickness 
increases rapidly with half spreading rate from 0.5 to 8 km, whereas 
subridge seawater penetration, inferred from the depth of faults 
deduced from seismic hypocenters, decreases as spreading rate in- 
creases 46 (Figs 3 and 4). We applied similar procedures to reconstruct 



paleo-oceans. Plate tectonic reconstructions, based on plate bound- 
aries and finite Euler vectors from refs. 47 and 48, were calculated for 
several Chrons (see Supplementary Information) since the Late 
Oxfordian (154.3 Ma). 

The subridge distribution of isotherms is related to mantle's 
potential temperature, spreading rate and vicinity to transform off- 
sets. The 150°C subridge isotherm rises with spreading rate: it lies 
< 0.5 km deep below mid-segment points at fast ridges versus 
> 4 km at ultra-slow ridges. At ridge-transform intersections it dee- 
pens depending on the age contrast, i.e., offset length and slip rate. 
One example is in the equatorial MAR where a mantle thermal 
minimum, enhanced by a long transform "cold edge effect", deter- 
mined nearly amagmatic accretion and an ultramafic seafloor along a 
~ 40 km long ridge segment 49,50 . 

Exposure of MORP at the seafloor, with the possibility of perido- 
tite/seawater reactions, is favored by vertical tectonics triggered by 
low-angle detachment faults and core complexes, common in slow- 
spreading ridges 40,51 " 53 , as well as by transtension/transpression along 
transforms due to small changes in ridge/transform geometry 38,54 . 
These processes may affect significant portions of the MAR with 
peridotites making up as much as 75% of the deeper rocks expo- 
sed 52,53,55 . In addition, significant stretches of slow ridges (i.e., Gakkel 
Ridge, SWIR, MAR) expose ultramafics along axis, away from trans- 
forms and core complexes 48,49,56,57 . Overall, up to 20% of the seafloor 
at modern mid ocean ridges may be floored by peridotites. 

These factors (Fig. 4) suggest that mantle-derived MORP can 
interact with seawater only at half spreading rates < 12 mm/yr (other 
than close to long transform intersections), i.e., in stretches of the 
mid Atlantic, Indian, American- Antarctic and Gakkel Ridges, but 
not along the fast-spreading East Pacific Rise (Fig. 2a), except in a 
few peculiar tectonic settings, i.e., Hess Deep 58 and the Garret 
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Figure 4 | Zone below mid ocean ridges where mantle-seawater reactions 
are possible. Basaltic crustal thickness increases with spreading rate, while 
depth of seawater penetration, inferred from earthquake hypocenters 
(yellow filled circles), decreases with spreading rate. The subridge zone 
with potential MORP-seawater reactions is indicated. Serpentinization, 
limited by the 500°C isotherm, is possible (other than close to long offset 
transforms) in ridges with half spreading rate < 12 mm/yr. MORP- 
seawater reactions at <150°C, where Mg-release is possible, are limited to 
ridges with half spreading rate < 8 mm/yr. Predicted crustal production 
and depth of the isotherms, obtained simulating two plates diverging along 
a mid ocean ridge, are indicated as red solid and dashed lines, respectively. 

Transform 59 . Adding up the slow-ultraslow-spreading ridge 
stretches, we estimate that approximately 9.2 - 10 u kg of mantle- 
derived MORP can interact yearly with seawater at T < 150°C and 
can potentially release Mg to the oceans. In places along slow mid 
ocean ridges, parts of layer 3 may be serpentinites instead of gabbros, 
both with a P-wave seismic velocity of ~ 6 km/s (Supplementary Fig. 
S3). Exposure of (sub-continental) mantle peridotites occurs also in 
magma-poor ocean-continent transition, as across the W-Iberian 
margin 60 . The zone of serpentinized peridotite along the W-Iberian 
margin is up to 6 km thick 60 , suggesting important alteration and 
serpentinization at depth. 

We estimated next the quantity of Mg contributed to the modern 
ocean by MORP-seawater reactions. If 100% of MORPs that can 
potentially react with seawater at < 150°C do actually react and lose 
at least 5%wtofMgO, seawater will gain 1.15" 10 12 mol/yr Mg. This 
is ~ 20% of the yearly Mg river input. If only 10% react, seawater will 
gain 1.15 • 10 11 mol/yr Mg, i.e., ~ 2% of the river input. Mg release by 
MORP-seawater reactions requires high water/rock ratios. Assum- 
ing W/R ratios from 10 to 10,000, the fluxes of H 2 0 required to 
mobilize 1.15- 10 12 mol/yr Mg (100% of MORPs) range from 
~ 10 13 to 10 16 kg/yr (i.e., close to or higher than the ridge high-T 
hydro thermal flux in basaltic crust). This is consistent with the total 
length of modern slow/ultraslow spreading ridges exceeding the 
length of fast spreading systems, and with larger volumes of rock 
being involved in low-T hydrothermal circulation. 

Moreover, ridge stretches where MORPs prevail imply a decrease 
of the ridge length where hydrothermal seawater reacts with basalt, 
lowering both the quantity of Mg extracted from seawater, and of Ca 



released to seawater (Fig. lb and lc). Although, low-T MORP-sea- 
water reactions at ridge flanks may be significant sources of Mg, we 
did not consider them, as they are difficult to quantify. We also did 
not consider contributions from non- volcanic passive margins, and 
supra-subduction zones. 

Scarcity of peridotite-seawater reactions in the Mesozoic Mega- 
ocean. Let us consider if and how the factors outlined above were 
different in the late Mesozoic, to see if they can explain the low Mg/ 
Ca ratio of Mesozoic seawater. Continents assemble in a single 
supercontinent and then gradually disperse again in ~ 500 million 
years Wilson cycles. A single supercontinent (Pangea) dominated the 
scene in the early Mesozoic 61 . This implies a single super-ocean 
(Panthalassa), a sort of super-Pacific, but hardly any intercon- 
tinental, Atlantic/Indian/ Arctic-type oceans. 

The super-ocean was at first probably almost surrounded by sub- 
duction boundaries and must have involved a fast-spreading "hot" 
mid ocean ridge system, generating a > 5 km thick basaltic crust, 
shallow 150°C and 500°C isotherms and shallow hydrothermal pen- 
etration, with little or no chance for mantle ultramafics to interact 
with seawater and to release Mg to the ocean. On the contrary, the 
Mesozoic super-ridge was probably the locus of intense high tem- 
perature hydrothermal circulation within the thick basaltic crust, 
with Mg being drained from, and Ca contributed to, seawater. 
Tethys, the only significant Cretaceous intercontinental ocean, 
hosted ridges that mostly spread faster than 20 mm/yr 62 , thus with 
no significant MORP-seawater reactions (Fig. 5). 

Seawater Mg/Ca temporal variations are paralleled by changes of 
8"B (ref. 24). 10 B is extracted from seawater preferentially to U B 
during hydration of mantle peridotites 63 . Thus, a U B increase in 
post-Mesozoic seawater, in parallel with the increasing importance 
of peridotite-seawater reactions, is consistent with our model. 

Mg isotopes may help in assessing the extent to which MORP- 
seawater reactions have contributed Mg to seawater. The 8 26 Mg of 
river input is today ~ — 1.09%o, different from that of seawater 
(-0.82%o) 64 . The 8 26 Mg of oceanic basalt is -0.36%o; given that no 
isotopic fractionation occurs during partial melting 65 , MORP should 
have a similar 8 26 Mg. Release of Mg from MORP to seawater, assum- 
ing no isotopic fractionation, should increase the 8 26 Mg of seawater 
relative to that of rivers, i.e., the difference between seawater and river 
water would become larger. In contrast, basalt-seawater interaction 
should not alter significantly seawater 8 26 Mg, assuming that Mg is 
entirely trapped by basalt or no isotopic fractionation during Mg- 
removal by hydrothermal circulation. We predict that the 8 26 Mg of 
Mesozoic seawater (little Mg input from MORP) was more negative 
than the 8 26 Mg of modern ocean (significant input from MORP). 
However, dolomite deposition, given that dolomite 8 26 Mg (- — 2%o) 
is more negative than that of rivers 64 , will also drive seawater 8 26 Mg 
towards less negative values, just as MORP-derived Mg would. 
Calculations on MORP-driven versus dolomite-driven changes of 
seawater 8 26 Mg suggest that modern ocean 8 26 Mg is affected signifi- 
cantly by MORP-derived Mg (see Supplementary Information). 

Weathering of continental rocks probably contributed a lower 
quantity of both Ca and Mg to the ocean during a "supercontinent" 
stage relative to a "dispersed continents" stage. A relatively low 
87 Sr/ s6 Sr ratio of Jurassic carbonates 66 supports this statement. 
However, a lower river input of Ca and Mg does not imply a different 
Mg/Ca ratio, given that the overall composition of the "dispersed 
continents" should not be very different from that of the "super- 
continent". 

Mg/Ca in post Mesozoic Oceans. Summing up, we surmise that the 
low Mg/Ca ratio of Mesozoic seawater was due not only to increased 
high-T hydrothermal circulation in basalt and to enhanced dolomite 
deposition, but also to the quasi-absence of seawater/mantle 
peridotite reactions in the Mesozoic. Significant Mg input to 
seawater due to hydration of mantle rocks must have started after 
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Figure 5 | Temporal variations of volume of mantle rocks that may react with seawater at T < 150°C along mid-ocean ridges since 150 Ma. 



the break up of the supercontinent and the gradual development of 
slow-spreading ridges in new intercontinental oceans. Based on plate 
tectonic reconstructions, a strong decrease in oceanic basaltic crust 
production and a sharp increase in MORP-seawater reactions 
occurred from the Santonian (anomaly M25, 83.5 Ma) to the late 
Paleocene (c25, 55.9 Ma) (Fig. 5 and Supplementary Information). 
Reconstruction at 83.5 Ma suggests that mid-ocean ridges half 
spreading rate was then generally faster than 15 mm/yr (Fig. 2b), 
with scarce possibility of low-T MORP-seawater reactions and of Mg 
release to seawater (Fig. 3). Different stretches of slow spreading ri- 
dges, with significant MORP exposures, became active at different 
times: Mg release to seawater must have increased accordingly. The 
Mid Atlantic Ridge/megatransform system, a major potential source 
of Mg, developed 80-100 Ma. The Gakkel Ridge-Lena Trough 
segments, with broad exposure of ultramafics 56,67 , developed about 
53 Ma and 10 Ma, respectively 68 . The SWIR from ~ 50° E to the 
Rodriguez triple junction, today a major MORP contributor, deve- 
loped from ~ 60 Ma to today 69 . The Andrew Bain megatransform in 
the SWIR, today a significant locus of ultramafic exposures, 
developed not earlier than about 50 Ma 70 . 

Discussion 

The cumulative increase of the quantity of MORP available to inter- 
act with seawater in the post-Mesozoic oceans (Figs. 1 and 5) paral- 
lels the increase of seawater Mg/Ca ratio documented by ref. 8. A 
>20% volume fraction of those MORPs potentially able to interact 
with seawater would help explain the Tertiary increase of Mg 
concentration in seawater (Fig. lc) inferred from halite fluid 
inclusions 3-5 . 

Model results suggest that three important processes have affected 
temporal variations of Mg/Ca during the last 150 Ma: (1) hydro- 
thermal circulation in basalts; (2) low- temperature MORP/seawater 
reactions; and (3) dolomite formation. Variations of oceanic crust 
production imply variations of Mg-removal and Ca-release by 
hydrothermal circulation in basalt. Given that the calculated oceanic 
crust production was higher during the Cretaceous than in the 
Cenozoic, we estimated a greater hydrothermal Mg-uptake and 
Ca-release in the Cretaceous than in the Cenozoic. We have shown 
that Mg- input to seawater due to Mg- release by MORP/seawater 
reactions was low (less than 4% of the river input) during the 
Cretaceous and high (up to 20% of the river input) during the 
Cenozoic. In addition, we have shown, from Mg isotope budget 



constraints, that Mg-capture by dolomite precipitation was higher 
during the Cretaceous (up to 23% of the river input) than during the 
Cenozoic (down to 14% of the river input). We conclude that the 
absence of Mesozoic MORP-seawater reactions, followed by their 
increasing importance, contributed significantly to the increase of 
seawater Mg/Ca ratio from the Mesozoic to the modern oceans. 

Methods 

Mg model. We assume a model where secular variations of seawater Mg-content are 
controlled by: (1) a constant influx from rivers; (2) variable Mg-release flux from 
peridotite- seawater reactions; (3) variable Mg-removal flux from high-temperature 
hydrothermal alteration of the basaltic crust as the result of seafloor spreading rate 
variations; and (4) a constant Mg residual-outflow including low-temperature off- 
axis hydrothermal interactions, deposition of dolomite, ion-exchange reactions with 
clays. Changes in the size of the oceanic Mg reservoir are thus calculated as: 

d ^=FZ g + ^ S p(t)-F™l(t)-F% (1) 

where: ivf? is the constant river influx, assumed of 5.6* 10 12 mol/yr; F^ p {t) is the 
influx due to MORP-seawater reactions. We let the MORP -derived Mg flux vary 
through time following the estimated volume of mantle rocks that can interact yearly 
with seawater at T < 150°C (Fig. 5). The estimated recent flux is of 1.15' 10 12 mol/yr, 
assuming that 100% of MORPs that can potentially react with seawater at T < 150°C 
do actually react and lose 5% wt of their MgO content (i.e., F^(0) = 0.05- p m -P m (0)/ 
[MgO mo i ar mass ], where p m — 3300 kg/m 3 is the density of mantle rocks and 
P m (0) = 2.7769" 10 s m 3 /yr is the current volume rate of MORP that interact with 
seawater at T < 150 U C); F^s is a constant unknown Mg-residual outflow; and 
Ff yB (t) = H hT (t) [Mg] sw (t) is the Mg-removal flux at time t by high-T hydrothermal 
circulation {Hf lT ) at ridge axis, that varies due to variations in the rate of seafloor 
production (Supplementary Tab. 55). [Mg] sw ( t) is the concentration of Mg in seawater 
at time t. Ht lT (0) is the modern mid ocean ridge high-T hydrothermal flux, assumed at 
5* 10 13 kg/yr. Assumed values of 3 and 5.2* 10 13 kg/yr (range of the estimated high-T 
hydrothermal flux) do not change our main results. 

Ca model. Secular variations of seawater Ca-content have been modelled assuming: a 
constant river influx F^ of 1.4* 10 13 mol/yr; a variable net influx due to Ca 

capture/release during low and high temperature hydrothermal circulation in mid 
ocean ridge basalts; a variable outflux F^orpW related* to MORP alteration at 
seafloor; and an unknown constant residual Ca-outflow F l r f s including: Ca-fixation 
due to carbonate accumulation (biogenic and inorganic), and anhydrite precipitation. 
Thus, changes through time of seawater Ca-content can be described by: 

^ =i£+J$W -F^MoAt) ~F% (2) 

We assume that secular variations of the net flux F^ B (t) due to Ca-release by 
hydrothermal circulation and Ca-capture by MORB alteration follow variations of 
oceanic crust production {Supplementary Tab. S5), with a modern MORB 
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hydro thermal- weathering net inflow F^ B (0) of 1.25" 10 12 mol/yr. Assumed values of 

0.8 ■ 10 12 and 1. 5 - 10 12 mol/yr (range of the estimated hydrothermal flux) do not 
change our main results. In addition, we assume variations through time of Ca- 
removal by MORP weathering F^ ORP (t) that scales linearly with variations in 
volume of mantle rocks that can interact yearly with seawater at T < 150 C {Fig. 5 and 
Supplementary Tab. S5). 

Numerical solutions of eqs. (1) and (2) were reached by finite difference 
approximation { Crank- Nicolson implicit scheme) using an integration time step of 
1 Ma and initial ( 1 50 Ma) Mg- and Ca-seawater concentrations ([Mg] = 30.5 mmol/ 
kg H 2 0 and [Ca] = 24 mmol/kg H 2 0), inferred from halite fluid inclusions 5 . The 
unknown residual fluxes F^f and F^f s were solved iteratively to fit the modern Mg- 
and Ca-seawater concentrations ([Mg] — 53 mmol/kg H 2 0 and [Ca] — 10.5 mmol/ 
kg H 2 0). Model details are presented in Supplementary Information. 
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